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ABSTRACT 

The magnification effects of clustered matter produce variations in the image sizes and number density 
of galaxies across the sky. This paper advocates the use of these effects in wide field surveys to map large- 
scale structure and the profiles of galaxy and cluster sized halos. The magnitude of the size variation 
as a function of angular scale is computed and the signal-to-noise is estimated for different survey 
parameters. Forthcoming surveys, especially well designed space-based imaging surveys, will have high 
signal-to-noise on scales of about 0.1 arcminute to several degrees. Thus the clustering of matter could 
be measured on spatial scales of about 50 Kpc to 100 Mpc. The signal-to-noise is dominated by sample 
variance rather than shot-noise due to the finite number density of galaxies, hence the accuracy of the 
measurements will be limited primarily by survey area, sampling strategy and possible systematics. 
Methods based on magnification are compared with the use of shape distortions and the contrasts and 
complementarities are discussed. Future work needed to plan survey strategy and interpret measurements 
based on magnification is outlined. 

Subject headings: cosmology: theory — gravitational lensing — large-scale structure of universe 



1. INTRODUCTION 

Gravitational lensing refers to the distortions in images 
of distant galaxies due to the deflection of light rays by 
mass concentrations. The distortion on a circular image 
can be decomposed into an amplification of the size of the 
image and an anisotropic stretching of its shape into an 
ellipse. The size amplification is called magnification and 
the anisotropic stretching is the shear. Gravitational lens- 
ing due to galaxy clusters and large-scale structure typi- 
cally leads to distortions of order 1-10% (e.g. Gunn 1967; 
Miralda-Escude 1991; Blandford et al 1991; Kaiser 1992; 
Bcrnardeau, van Waerbeke & Mellier 1997; Jain & Seljak 
1997; Kaiser 1998). In this regime of weak lensing the 
magnification /i is given by 

/i=((l- K ) 2 -|7| 2 ) _1 ^l + 2 K . (1) 

where k is the convergence and 7 the complex shear. So 
far observational studies of weak lensing have largely used 
the measured ellipticities to estimate the shear and thus 
the projected mass distribution. 

This paper makes the case for using effects of magnifi- 
cation in addition to the shear in mapping dark matter. 
Magnification leads to fluctuations in the sizes and num- 
ber densities (in a flux limited survey) of galaxies (e.g. 
Bartelmann & Narayan 1995; Broadhurst, Taylor & Pea- 
cock 1995; Schneider, King & Erben 200). In the context of 
galaxy clusters the change in number density has been used 
to constrain the mass distribution, but with less accuracy 
than the shape measurements (Taylor et al 1998). We ar- 
gue that for forthcoming blank field surveys the prospects 
are much better than for clusters to measure both effects 
of magnification, on sizes and number densities. 

(i) Magnification effects, unlike the shear, require con- 
trol fields to estimate the mean, unlensed size distribution. 
This had been a limitation for small, arcminute sized, clus- 
ter fields, but is automatically done in a blank field survey. 



(ii) The signal-to-noise (henceforth S/AT) due to shot 
noise is somewhat lower for magnification effects than the 
shear, and this has proven critical in cluster lensing, since 
a factor of 2 in S/Af is hard to make up. However field 
lensing surveys with areas larger than 10 square degrees 
are mostly in the regime where sample variance or sys- 
tematics dominate the errors. It is therefore feasible that 
from forthcoming imaging surveys with good control of 
systematics (photometric calibration for number density, 
resolution for sizes, and psf anisotropy for shear) all three 
lensing measurements can be made. Consistency checks on 
the different systematics can then be made, the S/Af on 
the measured dark matter clustering improved, and new 
information on halo properties can be extracted. 

(iv) Space based imaging surveys will make possible the 
measurements of sizes with an accuracy hard to achieve 
from the ground. Such surveys will become feasible over 
small areas with the Advanced Camera for Surveys on the 
HST, and over substantial fractions of the sky with a wide 
field imaging satellite telescope. 

The main goal of this paper is to propose that mea- 
surements of magnification effects, in particular the effect 
on galaxy sizes, be an integral part of the lensing agenda 
for forthcoming imaging surveys: wide area, multi-color 
ground based surveys like the CFHT Legacy Survey (see 
www.cfht.hawaii.edu/Science/CFHLS/), the proposed 
LSST (www.lss.org) and WFHRI (www.ifa.hawaii.edu/ 
"kaiser/wfhri) surveys, and especially a space based 
imaging survey as proposed for the SNAP satellite 
(http://snap.lbl.gov) which will have the key requirements 
of small psf and pixels ~ 0.1 arcsecond, photometric red- 
shifts, and survey area exceeding 100 square degrees (G. 
Bernstein, private communication). The formalism for 
computing statistical measures of magnification is pre- 
sented in Section 2. Section 3 provides estimates of the 
S/Af expected for measurements of size fluctuations. We 
conclude in Section 4. 
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2. STATISTICAL MEASURES OF MAGNIFICATION 

2.1. Fluctuations in the size distribution 

The lensing effect of an overdensity in the mass distri- 
bution is to increase the area of a given patch on the sky. 
The size of a given galaxy therefore increases. The area A 
and characteristic radius i?ocj4 1//2 ofa galaxy, for k <C 1, 
is then given by: 

A->A(1 + 2k): R^R(1 + k). (2) 

Thus the logarithm of the sizes is shifted linearly as 
log R — ► log R + k. 

Following the notation of Jain & Seljak (1997) 
we introduce the unperturbed metric ds 2 = 
o?{t) [-dr 2 + dx 2 + r 2 (d6 2 + sin 2 9d4> 2 )} , with r being 
conformal time, a the expansion factor normalized to 
unity today, \ the radial comoving distance and r(x) the 
comoving angular diameter distance. The convergence k 
is a weighted projection of the mass density along the line 
of sight. It can be expressed as 

*S) = ijfi™ r d x ffl ( x )^M, (3) 
^ Jo a 

where xh denotes the distance to the horizon. With W(x) 
denoting the radial distribution of galaxies in the sample, 
the radial weight function g(x) is given by 

g(x) = r(x) J™ r ^ { ~ f) x) w( x r )dx' ■ (4) 

The variance in the size fluctuations can be related to 
the variance in the smoothed convergence by consider- 
ing the mean size Rg in a circular aperture of angle 9. 
If the unlcnsed or intrinsic mean size in such an aper- 
ture is denoted Rg, then we can define the size shift 
5R = (Rg — R g )/R e . Using equations 2 and 3 the vari- 
ance of the shift in sizes along different lines of sight can 
be obtained in terms of the power spectrum of density 
fluctuations as 

{6R 2 }(9) = 36tt 2 fl 2 m / kdk / a" 2 (x) 
Jo Jo 

x Ps(k, X ) g 2 (x)U 2 [kr( X )0]dx, (5) 

where U(x) = 2J\(x)/x, with J\(x) being the Bessel func- 
tion of first order. 

2.2. Cross-correlations induced by magnification 

Magnification causes the observed area of a given patch 
of sky to increase, tending to dilute the number den- 
sity, but galaxies fainter than the limiting magnitude are 
brightened and may be included in the sample, thus in- 
creasing the number density. The net effect, known as 
magnification bias, can go either way depending on the 
slope s of the number counts of galaxies No(m) in a sam- 
ple with limiting magnitude m, s — d\ogN /dm. Mag- 
nification by amount fj, changes the number counts to 
(e.g. Broadhurst, Taylor & Peacock 1995), N'(m) = 
iVo(ra)/U 2 ' 5s-1 . In the weak lensing regime, this reduces to 
N'(m) = No(m) [1 + 5(s — 0.4)k] . Variations in the num- 
ber density which are correlated over some angular sepa- 
ration are produced due to the spatial correlations of the 



lensing dark matter. These correlations are difficult to de- 
tect since the galaxies have a strong auto-correlation func- 
tion due to their spatial clustering. However the cross- 
correlation of two galaxy samples with non-overlapping 
rcdshift distributions isolates the effect of magnification 
bias. 

The cross-correlation of a foreground-background 
galaxy sample can be obtained in the Limber approxima- 
tion (Moessner & Jain 1998): 

u x (9) = 3f> m (2.5s 2 - 1)4^ 2 dxW.ix) 9 -^ 

Jo a 

poo 

x / dkkP x ( X ,k) J [kr(x)9] , (6) 
Jo 

where the subscripts 1,2 denote the foreground and back- 
ground populations respectively and P x (Xi k) is the pro- 
jected galaxy-mass cross-power spectrum. oj x (9) can also 
be measured by the quasar-galaxy correlations that have 
been extensively discussed in the literature. The sizes of 
background galaxies discussed in the previous sub-section 
can be used as well; this would alter the equation above 
only in the numerical coefficient on the right-hand side. 

Note that u> x (9) is a measure of the galaxy-mass cross- 
correlation. It is the counterpart of galaxy-galaxy lens- 
ing, with the difference that the convergence is measured 
rather than the tangential shear. Hence it provides a more 
local measure of the galaxy-mass cross-correlation, which 
in the small scale regime probes the structure of galactic 
halos. For galaxy clusters one can measure size increases or 
number counts of background galaxies around individually 
clusters from high quality data, else they can be stacked 
like the galaxies. Large catalogs of clusters will soon be 
available for such measurements, and conversely, mass se- 
lected cluster catalogs may be obtainable from these mea- 
sures of the convergence. 

2.3. Corrections of high redshift supernovae magnitudes 

Recently Dalai et al (2002) have estimated the capabil- 
ity of shear maps to correct the lensing induced dispersion 
in the measured magnitudes of high redshift supernovae. 
The idea is the following. Magnification effects contribute 
to the measured scatter in the magnitudes of high redshift 
supernovae. Since the lensing contribution can equal or 
exceed the intrinsic dispersion for supernovae at z ^ 1, it 
is valuable to be able to measure the lensing effect along 
the lines of sight to individual supernovae from another 
tracer and thus correct the supernovae magnitudes. Dalai 
et al (2002) needed to assume that the convergence can 
be reconstructed on arcminute scales from ellipticity data, 
which may not prove to be feasible as they discuss. Mea- 
surements of size fluctuations however could directly map 
the convergence on arcminute scales around the line of 
sight to supernovae, allowing for a reduction in the scatter 
of supernova magnitudes. 

The main open questions are: How large is the shot 
noise effect for a given survey? How strongly is the 
smoothed convergence estimated from source galaxies cor- 
related with the value along the line of sight to individual 
supernovae? On arcminute scales ray tracing simulations 
can be used to estimate how much stronger this correlation 
is than assumed in the Gaussian limit taken by Dalai et al 
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(2002). Hence a quantitative study is merited to check if 
one can do better than the reduction of about 10% in the 
lensing dispersion reported by these authors. 

3. SIGNAL-TO-NOISE ESTIMATES 

We use as the primary observable the log of the linear 
size (such as half-light radius). To measure a lensing sig- 
nal, it is best to select galaxies on surface brightness which 
is conserved by lensing (Bartclmann & Narayan 1995). In 
the following we will assume that photometric redshifts 
are measured for the galaxies and that the surface bright- 
ness is not contaminated by atmospheric or instrumental 
effects. The number density of usable galaxies will depend 
on how conservative the catalog selection for a specific in- 
strument and survey will need to be. The S/N in the 
variance on smoothing scale 9 using galaxy sizes as an es- 
timator of the magnification is then given in terms of the 
standard deviation a% in the intrinsic solid angles and the 
number of galaxies N$ per circle of size 9. For a single 
field of size 9 the S/N is: 



1 -field 



a 2 M 



(7) 



where <J 2 (9) = (k 2 ) is the variance in the convergence with 
a top- hat smoothing of angular size 9, and <Ji{I) is the 
standard deviation of the size distribution in a given bin 
in physical surface brightness. This differs from the corre- 
sponding expression for ellipticity measurements because 
the denominator contains a different a 2 /Njj (the signal due 
to the variance of n or 7 in the numerator is equal in the 
weak lensing regime). From ground based data, the S/N 
for the ellipticity is larger than for the size since psf smear- 
ing directly affects the size estimate and affects the shape 
only at second order. From space based data however it 
is hard to say a priori whether the S/N from shape mea- 
sures is higher (at least by larger than a factor of two) for 
a given survey area. In any case, as shown below, on large 
scales sample variance dominates the statistical errors. 

For cosmological measurements, the data size of inter- 
est is a wide field survey from which the variance of k 
can be measured over a range of angular scales. Thus 
there are two angles, the first denoted Oo gives the size of 
the survey, and the second is the angular scale on which 
the variance is measured, which we will continue to de- 
note 9. Let Nf = Qq/9 2 be the number of patches of 
size 9 used in the measurement of k 2 . Thus the total 
number of galaxies is N t — NfN$. In the following we 
will assume that the Nf patches are uncorrelated. The 
contribution to the measured variance due to sample vari- 
ance and the intrinsic scatter in the size distribution is: 
N = [V2 a 2 K {6) + <J 2 /N g ] /y/WJ, where the first term is 
the sample variance contribution assuming a Gaussian dis- 
tribution, while the second is the contribution from the 
intrinsic scatter in the size distribution. 

The S/N for the measured variance is then given by 



[V2 <jl(9)+a 2 /N e ] 



(8) 



The above estimate ignores the effect of the kurtosis on 
the sample variance and thus underestimates the sample 



variance on small scales. The effect can be estimated by 
using the results of Takada & Jain (2002) who find that 
the kurtosis parameter defined as S4 — (k 4 )/<t® = 3x 10 4 
between 1' < 6 < 10' and falls off on larger scales (see their 
Figure 9). In the sample variance contribution, the rele- 
vant ratio is the standard definition of kurtosis in statis- 
tics, (k 4 )/<7 4 , so we need to find the angular scales on 
which this ratio is of order unity. Over the scales of in- 
terest, (T^~3x 10~ 4 where 9 is in arcminutes (Jain 
& Seljak 1997). Hence we obtain (k 4 )/ct 4 ~ lO/0('), a 
simple expression that is sufficiently accurate for our pur- 
pose. Thus for 9 ^ 10', the kurtosis term is important and 
could increase the sample variance by up to a factor of two. 
However, we will see below that the shot noise term dom- 
inates the sample variance term on scales smaller than a 
few arcminutes. Hence it is only over a small range in 
angle, and at worst by a factor of two, that we have un- 
derestimated the sample variance. Note that analogous 
expressions to equation 8 hold for the S/N from shape 
measurements (e.g. Jain & Seljak 1997; Schneider et al 
1998; Hu & Tegmark 1999). 

It is interesting to consider the relative contributions of 
sample variance and intrinsic scatter to the noise term 
in equation 8. Again using the approximate relation 
<r 2 cx 1/9, we see that the shot-noise term scales 1/9 3 , 
while the sample variance scales as 1/9 2 . Thus on small 
scales the shot-noise term dominates, while on scales larger 
than a few arcminutes (depending on the number density 
of galaxies) the sample variance term dominates. 

Figure 1 shows the predicted variance in size shifts and 
the S/N expected for different survey parameters. We as- 
sume a flat ACDM model with as — 0.9 and assume that 
photometric redshifts are available for a source redshift 
distribution of the form n(z) oc z 2 exp [— (z/zo) 1 ' 2 ] . Vary- 
ing zo changes the mean redshift of the distribution. The 
left panel shows the variance in the size shift and the two 
sources of noise: the intrinsic dispersion of galaxy sizes and 
sample variance. On scales of order 1' and smaller, the in- 
trinsic dispersion dominates, while on larger scales sample 
variance is the main source is noise. In the right panel the 
S/N achievable with a filled survey is shown. The middle 
solid curve assumes a number density of galaxies of 40 per 
square arcminute, a total area of 100 square degrees and 
an intrinsic dispersion Ui — 0.5 (Narayan & Bartelmann 
1995). This curve shows that high S/N measurements can 
be made on scales of order 0.1' — 100', which corresponds 
to spatial scales of about 50 Kpc to 50 Mpc. If the level 
of systematics is not a show-stopper, then one can extend 
the measurements to larger scales by sparse sampling. For 
given survey area, sparse sampling would increase the Nf 
term on large scales. Kaiser (1998) uses a power spectrum 
analysis to examine the best strategy for sparse sampling. 

The lower and upper solid curves in the right panel of 
Figure 1 show the effect of changing the survey area by a 
factor of ten — the curves shift up and down by the square 
root of the area. The dashed and dot-dashed curves show 
the effect on the S/N of lower galaxy number density and 
higher mean redshift of source galaxies, respectively. If 
the effective number density for which sizes can be mea- 
sured is decreased to 20 per square arcminute, then the 
shot noise term on small scales (9 < 5') lowers the S/N. 
For a higher redshift distribution of source galaxies, keep- 
ing other parameters constant, the signal is higher, so the 
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S/Af improves on small scales as shown by the dot-dashed 
curve. If neighboring fields are correlated, then the sam- 
ple variance estimate must be revised because the effec- 
tive number of independent fields of given angular size 9 is 
smaller. As discussed above, this would lower the S/Af for 

9 < 10'. It is clear from the range of the effective parame- 
ters in the S/Af explored here that even in a conservative 
scenario, a survey with area of order 100 square degrees 
will provide high S/Af measurements over several decades 
in length-scale. 

4. DISCUSSION 

What kind of survey would be suitable for measuring 
the magnification effects discussed in this paper? For the 
effect of magnification on galaxy sizes, a wide area space 
based multi-color imaging survey would be ideal. It is 
challenging for a ground based telescope to overcome the 
effect of psf smearing on the size distribution, unless one 
has the luxury of a large enough sample of galaxies with 
sizes larger than the psf. With appropriate multi-color 
imaging one can obtain photometric redshifts which can 
help reduce the scatter in measuring the size variance in- 
duced by lensing, allow one to check for intrinsic correla- 
tions in sizes and eliminate their contribution if needed. 
It also allows for the possibility of measuring the evolu- 
tion of matter clustering by binning the source galaxies in 
redshift (Jain & Seljak 1997; Hu 1999). With a psf of or- 
der 0.1 arcseconds and deep imaging, it is feasible to make 
size measurements on of order a million galaxies over a 

10 square degree area (based on the size vs. magnitude 
measurements in the Hubble Deep Field by Gardner & 
Satyapal 2000). This would give adequate S/Af to mea- 
sure the variance of the size distribution over a few bins 
in angle ranging from 1 to 10 arcminutes. 

With an area coverage of 100s of square degrees, which 
would probably be feasible only with a dedicated imag- 
ing satellite such as SNAP, one can measure the projected 
matter power spectrum to a precision of a few percent, 
measure higher order correlations, and ideally in combi- 
nation with shear information, get useful constraints on 
cosmological parameters. On the smallest scales, cross- 
correlation statistics would probe galaxy halos on scales 
of a few 10s of Kpc. By combining the magnification 
measurements with the shear, the density profile of ha- 
los can be measured far more accurately than with just 



galaxy-galaxy lensing, which probes only the integrated 
mass within radii. Further work is needed to quantify 
this, explore how small the scales that one can probe are, 
and check the validity of the approximation of equation 1 
on these scales. Magnification effects make possible other 
useful measures of the non-Gaussian lensing field that have 
proven difficult to obtain from shear data, such as the 
skewness of n which probes Q m (Bernardeau et al 1997) 
and peak statistics which probe the mass function of halos 
(Jain & van Waerbeke 2000). 

The cross-correlation effects of magnification on the 
number densities of galaxies, and of foreground galaxy po- 
sition with background galaxy sizes, are in principle easier 
to measure. This is because these statistics are first or- 
der in the lensing convergence whereas the size variance 
is of second order. The interpretation is more complex in 
that it involves the relation of a foreground galaxy pop- 
ulation with the mass. The main requirements for accu- 
rate measurements are photometric redshifts for a large 
sample of galaxies (to separate the foreground and back- 
ground galaxies), and high imaging quality as discussed 
above. For deep imaging data that has a redshift distri- 
bution peaked at z 1, an adequate dataset would en- 
compass 10 square degrees, while an ideal dataset would 
cover more than a 100 square degrees. The southern strip 
of SDSS fulfills the requirements outlined above, as do 
other smaller imaging surveys that are in progress or be- 
ing planned. 

It is hoped that the discussion and results presented here 
motivate the integration of magnification measurements as 
part of the scientific agenda of wide area imaging surveys. 
The precise requirements for a given survey that will en- 
able useful magnification measurements to be made need 
careful consideration. At the same time work is needed on 
survey strategy, techniques for combining magnification 
and shear information, and appropriate statistical mea- 
sures that can be extracted from the data. 
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Fig. 1. — (a) The left panel shows the variance in fluctuations in the sizes of galaxy images vs. angle for a 100 square degree survey. The 
solid curves shows the predicted variance induced by lensing for a ACDM model for mean source redshift z s = 1 (lower curve) and z s = 2 
(upper curve). The dashed curves shows the contribution due to intrinsic size dispersion for galaxy number density n g = 40 per square 
arcminute (lower) and n g = 20 (upper). The dotted curve shows the sample variance contribution for z s = 1 (lower) and z s = 2 (upper), 
(b) The right panel shows the S/Af expected for a 10 square degree survey (lower solid curve), a 100 square degree survey (thick solid curve), 
and a 1000 square degree survey (upper solid curve) with n g = 40 and z s = 1. The dashed and dot-dashed curves show the decrease in S/Af 
for the middle curve if only half the number density of galaxies is available, and the increase if z s = 2, respectively. Sparse sampling would 
enhance the S/Af for large 8 compared to the results shown here which assume a filled survey. 



